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ABSTRACT 

There are currently two well-accepted models that explain how pulsars exhibit glitches, sudden changes 
in their regular rotational spin-down. According to the starquake model, the glitch healing parameter, 
Q, which is measurable in some cases from pulsar timing, should be equal to the ratio of the moment of 
inertia of the superfluid core of a neutron star (NS) to its total moment of inertia. Measured values of the 
healing parameter from pulsar glitches can therefore be used in combination with realistic NS structure 
models as one test of the feasibility of the starquake model as a glitch mechanism. We have constructed 
NS models using seven representative equations of state of superdense matter to test whether starquakes 
can account for glitches observed in the Crab and Vela pulsars, for which the most extensive and accurate 
glitch data are available. We also present a compilation of all measured values of Q for Crab and Vela 
glitches to date which have been separately published in the literature. We have computed the fractional 
core moment of inertia for stellar models covering a range of NS masses and find that for stable NSs in 
the realistic mass range 1.4 ± 0.2 Mq, the fraction is greater than 0.55 in all cases. This range is not 
consistent with the observational restriction Q < 0.2 for Vela if starquakes are the cause of its glitches. 
This confirms results of previous studies of the Vela pulsar which have suggested that starquakes are not 
a feasible mechanism for Vela glitches. The much larger values of Q observed for Crab glitches (Q Z 0.7) 
are consistent with the starquake model predictions and support previous conclusions that starquakes 
can be the cause of Crab glitches. 

Subject headings: pulsars - individual: Crab (B0531+21), Vela (B0833— 45) 



1. INTRODUCTION 

Soon after the discovery of radio pulsars (Hewish et al. 

1968) , they were identified as rotating, highly magnetized 
neutron stars (NSs) (Gold 1968, 1969; Ostriker & Gunn 

1969) . It was observationally established that pulsar peri- 
ods slowly and monotonically increase with time as a re- 
sult of magnetic braking and rotational energy loss. Early 
in 1969, it was noticed that the Vela pulsar suddenly in- 
creased its angular velocity (Reichley & Downs 1969; Rad- 
hakrishnan & Manchester 1969). These distinct and sud- 
den increases in rotational frequency, known as glitches, 
were later regularly observed in the Vela and Crab pul- 
sars and, more infrequently, in other pulsars (e.g., Shemar 
& Lyne 1996; Lyne, Shemar, & Smith 2000). The Crab 
(PSR B0531+21) and Vela (PSR B0833-45) pulsars are 
the two best-studied glitching pulsars since they have each 
been observed to glitch a number of times and are bright 
and easy to monitor. These pulsars have the most ex- 
tensive and accurate glitch data published and currently 
provide the best test of the physical mechanisms by which 
pulsar glitches occur. 

There are currently two competing models that are well- 
accepted to explain how pulsars glitch. In simplest terms, 
both models treat the NS as a two-component body, with 
a superfluid interior core surrounded by a rigid external 
crust (e.g., Ruderman 1972). These two components are 
dynamically weakly coupled through the magnetic field. 



In the starquake glitch model (Ruderman 1976; Baym & 
Pines 1971; Alpar et al. 1996), a slight equatorial oblate- 
ness in the crust can be formed if the NS is born rapidly 
spinning. As the pulsar slows down via magnetic braking, 
the deformation cannot be supported through centrifugal 
pressure from the core, and the crust can suddenly crack 
under gravity. The subsequent reduction in oblateness re- 
duces the moment of inertia, resulting in a sudden increase 
in the rotational frequency, seen as a glitch. 

Alternatively, in the vortex unpinning model of glitches 
(Anderson & Itoh 1975; Alpar et al. 1984, 1993), angu- 
lar momentum is stored in vortices of superfluid which 
migrate outward from the core as the star slows down. 
These vortices can become pinned to nuclei in the transi- 
tion region between the core and outer crust, thereby pre- 
venting angular momentum transfer to the crust. A differ- 
ential rotation develops between the core and crust until 
a catastrophic unpinning of the vortices occurs. Angular 
momentum is suddenly transferred to the crust, spinning 
it up. Since the crust is tightly coupled to the external 
magnetic field, this sudden spin-up in the rotation period 
is observed as a glitch. 

Within the starquake model, the time evolution of a pul- 
sar's rotational frequency after the occurrence of a glitch 
(Wang et al. 2000) can be fit by the equation: 



n(i) = Q (t) + Af7 t=0 1 - Q(l - e 
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Here £lo(t) represents an extrapolation of the pre- glitch 
frequency evolution, Af2 t= o is the magnitude of the change 
in the rotational frequency at the time of the glitch (t = 0), 
t is a characteristic exponential healing or recovery time 
after the glitch, and Q is the fraction of the initial fre- 
quency change that is eventually recovered (known as the 
healing parameter) . The last term accounts for any perma- 
nent (not recovered) change AQ p in the frequency deriva- 
tive after the glitch. The parameter Q in Equation 1 can 
be determined observationally (Shapiro & Teukolsky 1983) 
from the jump in and its first and second time deriva- 
tives: 



Q = 



(2) 



t=o 



In the starquake model, Q is related to the moments of 
inertia of the components of the star according to (Pines, 
Shaham, & Ruderman 1974): 



Q 



/total 



(3) 



where I core and /total are respectively the moment of in- 
ertia of the superfluid core and the moment of inertia of the 
entire star. We have constructed a series of models of the 
NS interior using representative high-density equations of 
state (EOSs) in order to test whether the calculated mo- 
ments of inertia in the models satisfy the prediction of 
the starquake model for Crab and Vela glitches. We have 
compared the model results with the measured values of Q 
from Crab and Vela glitches to see if they are consistent. 

2. EQUATIONS OF STATE 

We have constructed NS models using seven representa- 
tive EOSs of superdense matter from which parameterized 
model stars could be produced. For densities p < 1.6 x 10 14 
g cm~ 3 (crustal densities), we have used three EOSs for 
three separate density regimes (see, e.g., Shapiro & Teukol- 
sky 1983). These three EOSs are described by Feynman, 
Metropolis, & Teller (1949), Baym, Pethick, & Sutherland 
(1971), and Baym, Bcthc, & Pethick (1971). The seven 
high-density EOSs, which dominate the macroscopic char- 
acteristics of the model stars, are briefly described below. 
Further details about each EOS can be found in the sources 
referenced. 

2.1. The BJW Equation of State 

For densities above the nuclear density, a combination of 
EOSs presented by Canuto (1975) was used for the mod- 
eling. For densities up to ~ 5 x 10 15 g cm -3 , the EOS 
described by Bethe & Johnson (1974) was employed. This 
EOS is an improvement upon the work of Reid (1968) 
which includes a repulsive nucleonic core arising from the 
exchange of vector mesons in a hyperonic liquid. At higher 
densities, the EOS of Walecka (1974) was used, although 
the NS models produced in this density regime were be- 
yond the stability limit. In this EOS, nucleons interact 
attractively via exchange of scalar mesons and repulsively 
via exchange of more massive vector mesons. Summaries 
of these EOSs are also presented by Shapiro & Teukolsky 
(1983). We call this EOS BJW. 



2.2. The FPS Equation of State 

Lorenz, Ravenhall, & Pethick (1993) refer to an EOS 
calculated by Friedman & Pandharipande (1981) for high- 
density neutron and nuclear matter. This EOS employs 
the microscopic V14 two-body potential, a three nucleon 
interaction (TNI) potential, and uses hypernetted chain 
techniques. This EOS was modified by fitting the mi- 
croscopic interaction of Friedman & Pandharipande to a 
Skyrme-like energy density function (Skyrme 1959). The 
essential feature of the Skyrme model is a two-body inter- 
action that has the spatial character of a two-body delta 
function plus derivatives. We refer to this EOS as FPS. 

2.3. G WM Equation of State 

This EOS is based on a variant of the theory of nu- 
clear field coupling (Zimanyi & Moszkowski 1990) in which 
the scalar field is coupled to the derivative of the nucleon 
field. A modification proposed by Glendcnning, Weber, 
& Moszkowski (1992), called the hybrid derivative cou- 
pling model, replaces the purely derivative coupling of the 
scalar field to baryons and vector mesons with a Yukawa 
point and derivative coupling to baryons and both vector 
fields. The coupling model is consistent with the exper- 
imentally inferred binding energy of lambda hyperons in 
nuclear matter. The resulting EOS, called GWM here, 
maintains equilibrium between all baryons to convergence 
and leptons. The existence of hyperons significantly soft- 
ens this EOS. 

2.4. HKP Equation of State 

This EOS, proposed by Hacnscl, Kutschera, & Proszyn- 
ski (1981), supposes that a minimal condition for the prop- 
erties of high-density matter is that it must produce the 
observed values for the saturation properties of nuclear 
matter. The relativistic mean field theory of Serot (1979) 
is assumed in which nucleon interactions are governed by 
the exchange of neutral scalar and neutral vector mesons, 
pions, and rho mesons. This theory is an extension of the 
relativistic mean field theory of Walecka (1974). The au- 
thors assume a conservative range of possible nuclear sat- 
uration densities, the mean of which is 2.8 x 10 14 g cm~ 3 . 
The EOS used here, called HKP, is based on this nuclear 
saturation value. 

2.5. WFF Equation of State 

Wiringa, Fiks, & Fabrocini (1988) describe a mi- 
croscopic EOS of dense nuclear matter constrained by 
nuclcon-nucleon scattering data. The interaction in- 
cludes a two-nucleon Urbana V14 (UV14) potential and a 
TNI three-body potential term (Lagaris & Pandharipande 
1981a, b). The three-body potential includes a repulsive 
term, the primary effect of which is a reduction in the 
intermediate-range attraction of the two-nucleon poten- 
tial, and an attractive term, which becomes negligible at 
high densities. We refer to this EOS as WFF. 

2.6. APR Equations of State 

Akmal, Pandharipande, & Ravenhall (1998) describe a 
set of realistic EOSs based on the Argonne v\s (A18) two- 
nucleon interaction (Wiringa, Stoks, & Schiavilla 1995), 
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calculated using variational chain summation methods. 
This is supplemented with a relativistic boost correction 
term (Forest, Pandharipande, & Friar 1995) and a TNI 
term based on the Urbana IX model (Pudliner et al. 1995). 
The two cases of pure neutron matter (PNM) and symmet- 
ric nuclear matter (SNM) are separately considered. The 
latter is composed of equal numbers of neutrons and pro- 
tons in beta equilibrium. The two EOSs used to construct 
the models presented here were derived from a parabolic 
fit to and extrapolation of the tabulated nucleon density 
and energy values given for PNM and SNM by Akmal, 
Pandharipande, & Ravenhall (1998). We refer to the two 
EOSs based on PNM and SNM as APR(p) and APR(s), 
respectively. 

3. THE MODELING PROCEDURE 

To construct parameterized NS models for each EOS, 
we first chose a central density p c = p(0) for each stel- 
lar model. The EOS and the relativistic expressions for 
stellar structure (Oppenheimer & Volkoff 1939) were used 
to calculate the pressure gradient and incremental mass 
contained in concentric shells as the model iterated out- 
ward from the center of the star. The relativistic structure 
equations are: 



dm{r) 
dr 



= AiTr 2 p(r) 



(4) 



dPjr) 
dr 



p{r) 



P(r)\ G [m(r) + Airr 3 P (r) / c 2 ] 
r 2 [l-2Gm(r)/rc 2 ] 



(5) 

where m(r) represents the mass of the star internal to a 
distance r from the center of the star. When the pressure 
reached zero, P(R) = 0, the edge of the star was reached 
(r = R) and the relevant macroscopic parameters could be 
read off. 

To calculate the relativistic moment of inertia I(r) in- 
terior to a radius r, expressions given by Arnett & Bow- 
ers (1977) were used which account for the Lense-Thirring 
frame-dragging effect (e.g., Glendenning 2000) (see discus- 
sion below): 



dl(r) _ 8tt [p(r)+P(r)/c 



21 A - 

z 'r e 



-0(r) 



dr 



ft 



(6) 



3 [l-2Gm(r)/rc 2 ] 1 /2 

Here il is the observed stellar rotational frequency and 
u>(r) is the angular velocity of the star at a distance r rel- 
ative to the angular velocity of the rotating local inertial 
frame, which is dragged at angular velocity oj(r): 



uo{r) = ft — uj(r). 



(7) 



The term e~ < ^ r - ) in Equation 6 translates u>(r) as mea- 
sured from infinity to u)(r) as measured in the local iner- 
tial frame at r. The following auxiliary relations were used 
to determine values of <p(r) and u{r) throughout the star 
(r < R) in the moment of inertia calculation: 



d<P(r) 
dr 



1 



p{r)c 2 + P{r) 



dP{r 
dr 



1 2Gm(r) \ ^ 



(8) 
(9) 



d_ 

dr 



iz ( r4 ?( r ) 



du>(r) 
dr 



4r 3 ^ 
dr 



Cu{r) = 0. (10) 



An arbitrary value of w(0) was chosen and the following 
boundary conditions imposed in the models: 



tfc)(0) 
dr 







u>{R) = Ct — 



RdCo(R) 
3" dr ' 



(11) 



(12) 



This latter condition (in which a value for fi is deter- 
mined at the end of the model iteration) was used to scale 
the arbitrarily chosen u>(0) in such a way as to yield the 
proper observed O in a subsequent modeling pass. This 
subsequent pass produced the correct u)(r) for all r < R. 
The condition of finite </>(0) was also imposed. 

Milne's centered algorithm (Harrison et al. 1965) was 
used for the iteration. Since there is a transition region 
between the superfluid core and the rigid crust which con- 
sists of a mixture of nuclei and free baryons, there is 
no clear division where the core-crust transition occurs. 
We chose a transition density near the nuclear density 
(/^transition ~ 2.4 x 10 14 g cm~ 3 ) in order to separate the 
two components in our models. From this we obtained 
^coro = ^transition) and J tota] = I (R) for each model NS 
and could compare the ratio of these to the measured val- 
ues of Q from glitch observations. In previous work, some 
authors (e.g., Datta & Alpar 1993; Shapiro & Teukolsky 
1983) have assumed a slightly lower core-crust transition 
density than the one used here. This lower transition den- 
sity reduces the size of the crust, thereby increasing the 
fractional moment of inertia of the core. Although slight 
differences in the choice of the transition density do not af- 
fect our model results, a significant reduction would lead 
to a more stringent test of the starquake model (i.e., it 
would further limit the range of observed Q that would be 
consistent with the predictions of the starquake model). 
Thus, our choice of /ftransition is a conservative one for this 
purpose. 

For angular rotation speeds comparable to those of the 
Crab and Vela pulsars (ft ~ 190 rad s" 1 and 70 rad s _1 , 
respectively) , the macroscopic star parameters are affected 
by first-order rotation effects, but not significantly so by 
second-order effects. 

The Lense-Thirring frame dragging effect is a first-order 
effect which scales with rotational speed as (Q/Q c ) (Arnett 
& Bowers 1977), where n c = (GM/R 3 ) 1 / 2 , the critical 
value above which centrifugal mass shedding would occur 
for the rotating star. For realistic NS models, fl c ~ 13000 
rad s _1 . This first-order effect is therefore ~ 1-2% for the 
Crab and Vela pulsars and is taken into consideration in 
the calculation of the relativistic moment of inertia in the 
models. 

NS mass shifts and deformations in sphericity from rota- 
tion (which would complicate the simple spherical stellar 
modeling presented here) are second order effects, scaling 
as (£l/£l c ) 2 (Arnett & Bowers 1977). This term would be 
significant for millisecond pulsars, but for the Crab pulsar, 
the second-order effect is ~ 10~ 4 , and it is even lower for 
Vela. Thus, the second-order effect can be safely ignored 
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when considering the macroscopic parameters of our NS 
models. 



4. RESULTS 

We have produced a set of parameterized NS models us- 
ing the seven high-density EOSs described above. In order 
to decide which models in each series were appropriate to 
use in the comparison with the observed values of Q, we 
used a realistic NS mass range as a reasonable constraint. 
Models with a NS mass outside of the range were not used 
in the comparison. 

Thorsett & Chakrabarty (1999) have presented a sta- 
tistical study of NS masses for the known binary radio 
pulsar population for which useful mass constraints could 
be derived. They find that, statistically, the distribution of 
pulsar masses falls in the range M = 1.35 ±0.04 M©. This 
is consistent with the initial NS mass function from simu- 
lations of supernovae, which predict that NSs are formed 
with M> 1.2 M© (Timmes, Woosley, & Weaver 1996). 
For our comparison and analysis, we consider only the 
mass range M = 1.4 ± 0.2 M©, which we assume to be 
a conservative estimate of the realistic NS mass range. 

Figure 1 shows NS mass M = m(R) as a function of cen- 
tral density p c = p(0) for our models constructed using the 
seven different EOSs described above. NS branches with 
dM/dp c > satisfy a necessary stability criterion (e.g., 
Arnett & Bowers 1977) and are indicated with solid bold 
lines. For all EOSs used, there are stable model stars in 
the mass range 1.4 ± 0.2 M©. Figure 2 shows NS mass M 
as function of radius R for the seven EOSs used. Stable NS 
branches are again indicated with solid bold lines. None 
of the models in the stable branches violates the required 
stability condition M/R < 4c 2 /9G for static relativistic 
stars (Glendenning 2000). 

Figure 3 shows the ratio of the moment of inertia of 
the core, I CO ro, to the total moment of inertia, /total, as a 
function of NS mass M. For stable models in the mass 
range of interest (1.4 ± 0.2 M©), this ratio is confined to 
^coro/^totai ~ 0.55. Starquakes cannot be responsible for 
producing glitches in a pulsar with measured values of Q 
consistently outside this range unless the pulsar mass is 
significantly smaller than expected. 

Glitch measurements (including measurements of Q) for 
the Crab and Vela pulsars have been published in various 
places in the literature (see Table 1 for references). How- 
ever, no comprehensive current listing of all measured val- 
ues of Q for Crab and Vela glitches exists. Table 1 is a 
compilation of all measured values of Q published to date 
for glitches from the Crab and Vela pulsars. Some are 
different measurements of the same glitch. This complete 
set of measurements can be used to derive a range of ob- 
served Q for each pulsar. From the 21 measurements of 
Q for Crab glitches, a weighted mean of the values yields 
Q = 0.72 ± 0.05 1 . This is only slightly smaller than the 
unweighted mean of Q = 0.83. A range of Q > 0.7 en- 
compasses the observed distribution for the Crab pulsar. 
A weighted mean of the 11 measurements of Q for Vela 
glitches yields 0.12 ±0.07, while an unweighted mean gives 
Q = 0.18. All estimates for Vela agree that Q is small, with 
a likely range Q < 0.2. 

1 For measurements of Q without a quoted uncertainty, an uncertainty 



5. DISCUSSION 

The glitch behavior of the Crab (Alpar et al. 1994, 1996) 
and Vela (Alpar et al. 1993; Chau et al. 1993) pulsars has 
been extensively studied previously. Below we discuss and 
compare our model results with the observed range of Q 
in the context of this previous work. 

5.1. Vela 

Alpar et al. (1993) and Chau et al. (1993) have found 
that their model predictions of vortex unpinning in Vela 
are completely consistent with the observed Vela glitch 
characteristics. They deduce that the fractional crustal 
moment of inertia must be greater than 2.6% for Vela 
(Chau et al. 1993). Datta & Alpar (1993) have used a 
corresponding lower limit of 3.4% for Vela in an attempt 
to rule out soft EOSs. Link, Epstein, & Lattimer (1999) 
employ a lower limit of 1.4% for a similar purpose. These 
estimates imply an upper limit to the fractional core mo- 
ment of inertia of ~ 0.98, which tends to disfavor APR(s) 
models with mass M £ 1.5 M©, but cannot independently 
constrain the Vela mass since this condition is satisfied for 
a wide range of model masses (see Figure 3). The conclu- 
sions made by Alpar et al. (1993) and Chau et al. (1993) 
that the Vela mass is probably less than 1.4 M© (and is 
more likely closer to 1.2 M©) cannot be confirmed or ruled 
out. Lyne (1992) explains one reason why starquakes can- 
not be the glitch mechanism for Vela: the required oblate- 
ness is not sustainable given the large size of Vela glitches 
(Afi/fi ~ 10~ 6 ). After — 100 years, the oblateness would 
reach zero, and therefore glitches could only have been 
sustained for ~ 1% of the current Vela age. The time in- 
terval between Vela glitches is also inconsistent with the 
much longer intervals predicted by the starquake model 
(Alpar et al. 1996) unless a solid core model is invoked for 
Vela (e.g., Canuto & Chitre 1973; Pines, Shaham, & Ru- 
derman 1974). Starquakes would also produce a thermal 
energy dissipation during the large Vela glitches which is 
expected to be observable as a change in X-ray luminosity 
soon after the glitch occurs. X-ray observations of Vela 
show no such signal down to limits of less than a few per- 
cent in the fractional change in flux (Seward et al. 2000; 
Helfand, Gotthelf, & Halpern 2001). The vortex unpin- 
ning model does not have this energy dissipation problem 
for Vela (Alpar 1995). The low values of Q measured for 
Vela glitches are additional evidence that the predictions of 
the starquake model do not match observations; otherwise 
the implied Vela mass would be too low (M < 0.5 M© for 
Q $ 0.2 in our models). Our model results confirm previ- 
ous conclusions that starquakes cannot be a feasible glitch 
mechanism for the Vela pulsar. 

5.2. Crab 

Alpar et al. (1994, 1996) have studied the behavior of 
Crab glitches and have determined that the starquake 
model is at least partially responsible for them. The mag- 
nitudes of Crab glitches, combined with the observed glitch 
rate, support this notion. Lyne (1992) indicates that the 
small size of Crab glitches (Afl/fl ~ 10~ 8 ) allows for small 
changes in oblateness which would not significantly de- 
plete the oblateness over the current lifetime of the Crab. 

of 0.1 was simply assumed in the calculation of the weighted mean. 
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Link, Franco, & Epstein (1998) and Franco, Link, & Ep- 
stein (2000) also show that the permanent post-glitch off- 
sets in the period derivative seen for the Crab pulsar can 
be accounted for by starquakes: the net torque on the 
star increases through shearing effects. Our model results 
are also consistent with the starquake interpretation for 
Crab glitches. The fractional moment of inertia values in 
our models are comparable to observed Q values for Crab 
glitches (Q > 0.7) if M > 0.15 M , which seems likely. Al- 
par et al. (1994) set a lower limit of 0.2% for the fractional 
crustal moment of inertia. This is not constraining in our 
models since the corresponding fractional core moment of 
inertia range icore/^totai ~ 0.998 is easily satisfied. Our re- 
sults support previous suggestions that starquakes could 
be responsible for Crab glitches. 

6. CONCLUSIONS 

Using parameterized NS models produced from seven 
representative EOSs of superdense matter, we find that 
the fractional moment of inertia of the core component of 
the model stars is /core/^totai ~ 0.55 for all stable configu- 
rations in the assumed realistic NS mass range 1.4±0.2 M 
(Figure 3). This ratio, which is predicted by the starquake 
model to equal the glitch healing parameter Q, is not con- 
sistent with the observed range Q < 0.2 for Vela glitches 



(see Table 1), unless the Vela pulsar mass is unrealistically 
small (M < 0.5 M@). This confirms results from previous 
studies of the Vela pulsar which indicate that starquakes 
are not the cause of Vela glitches (Lyne 1992; Alpar et 
al. 1993; Chau et al. 1993; Alpar 1995). The much larger 
values of Q > 0.7 seen for Crab glitches (Table 1) are con- 
sistent with the moment of inertia values of our models 
for realistic masses, as predicted by the starquake model. 
These results support previous conclusions from the anal- 
ysis of Crab glitch behavior (Alpar et al. 1994, 1996; Link, 
Franco, & Epstein 1998; Franco, Link, & Epstein 2000) in 
which starquakes have been proposed as the Crab glitch 
mechanism. Repeated and accurate measurements of Q for 
other glitching pulsars in the future (if obtainable) could 
be used as a simple test of the starquake glitch model and 
may help resolve whether differences between Crab-like 
and Vela-like glitches can be understood on evolutionary 
grounds (e.g., Alpar 1995). 

We thank the referee for helpful comments to improve 
the manuscript and for suggesting the inclusion of addi- 
tional equations of state in the modeling. MD was par- 
tially supported by a grant of the Polish Committee for 
Scientific Research. 
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NS Central Density (g cm 3 ) 



Fig. 1. — Mass vs. central density for NS models constructed from seven high-density EOSs. The seven parameter curves correspond to 
these seven EOSs in order of increasing maximum mass: GWM, BJW, FPS, WFF, APR(p), APR(s), and HKP. Branches with positive slope 
(a necessary condition for stability) are indicated with bold lines. The dashed horizontal lines correspond to the mass range 1.4 ± 0.2 Mq, 
which is assumed to be a realistic NS mass range. 
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Fig. 2. — Mass vs. radius for NS models constructed from seven high-density EOSs. Stable branches are indicated with bold lines. The 
dashed horizontal lines correspond to the mass range 1.4 ± 0.2 Mq, which is assumed to be a realistic NS mass range. 




NS Mass (solar masses) 

Fig. 3. — Fractional core moment of inertia vs. mass for NS models constructed from seven high-density EOSs. The assumed realistic NS 
mass range 1.4 ± 0.2 Mq is indicated by the dashed vertical lines. Stable NS configurations (in bold) in the mass range 1.4 ± 0.2 Mq have 
corresponding /core /-/total ~ 0.55 in all cases. This ratio is consistent with the large values of the glitch healing parameter Q predicted by the 
starquakc model for Crab glitches (Q }Z 0.7), supporting the starquakc glitch interpretation. The ratio is inconsistent with the much lower 
values of Q seen for Vela (Q $ 0.2), indicating that starquakes do not account for Vela glitches. 
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Measured Values of Q for Crab and Vela Glitches. 
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